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Abstract 

AKARI/Infrared Camera (IRC) Point Source Catalog provides a large amount 
of flux data at S9W (9 /im) and L18W (18 /im) bands. With the goal of con- 
structing Star-Formation Rate(SFR) calculations using IRC data, we analyzed an IR 
selected GALEX-SDSS-2MASS-AKARI(IRC/Far-Infrared Surveyor) sample of 153 
nearby galaxies. The far-infrared fluxes were obtained from AKARI diffuse maps to 
correct the underestimation for extended sources raised by the point-spread function 
photometry. SFRs of these galaxies were derived by the spectral energy distribution 
fitting program CIGALE. In spite of complicated features contained in these bands, 
both the SOW and Li^VF emission correlate with the SFR of galaxies. The SFR 
calibrations using S9W and L18W are presented for the first time. These calibrations 
agree well with previous works based on Spitzer data within the scatters, and should 
be applicable to dust-rich galaxies. 
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1. Introduction 



The star formation activity is fundamental to studies on the formation and evolution of 
galaxies. Numerous efforts have been made to find reliable and convenient SFR indicators (e.g. 
Kennicutt, 1998; Hirashita et al., 2003; Bell, 2003; Hopkins et al, 2003; Calzetti et al., 2007). 
Among most frequently used indicators, the ultraviolet (UV) and the optical recombination lines 
(e.g. Ha, Paa) give direct measures of light from young stars. However, UV and Ha emissions 
are strongly affected by dust extinction (Kennicutt, 1998), and Paa as well as other optical 
recombination lines were also shown to underestimate SFRs for galaxies with high luminosity 
(Rieke et al., 2009). On the other hand, since dust absorbs UV/optical light and re-emits 
the bulk of the energy into far-infrared (FIR) band (25 /im ~ 350/im), FIR emission could 
efficiently trace SFRs for dusty galaxies. However, FIR emission is unable to trace the dust- 
unobscured radiation and includes part of radiation from old stellar populations. Therefore the 
energy balance method combining the FIR and UV derived SFRs were used to complement the 
emission from young stars not traced by FIR (e.g. Buat & Xu, 1996; Meurer et al., 1999; Gordon 
et al., 2000; Buat et al., 2005). Although this method could trace SFRs with considerable 
accuracy, it is difficult to obtain the total dust emission especially for high redshift objects. 

The MIR monochromatic fluxes were also investigated as SFR indicators. The MIR 
emission is contributed by several components, including the polycyclic aromatic hydrocarbon 
(PAH) features (prominent at 6.2, 7.7, 8.6, 11.3, 12.7 and 17 /im), the continuum by the 
stochastic heating of very small grains, silicate absorption at 9.7 and 18 /xm, molecular hydrogen 
lines and fine-structure lines (Leger & Puget, 1984; Leger et al., 1989; Desert et al., 1990; Draine 
& Li, 2007; Smith et al, 2007; Treyer et al., 2010). The MIR-SFR relation was intensively 
studied using Spitzer IRAC 8 /im and MIPS 24 /im photometry data and spectral data (e.g. 
Forster Schreiber et al., 2004; Calzetti et al., 2005, 2007; Kennicutt et al., 2009; Rieke et al, 
2009; Treyer et al., 2010). Nevertheless, there is still a debate about the reliability of MIR 
indicators because of the complicated features it contains. Attempts to combine the optical 
and IR indicators show that this combination could trace the SFR effectively, and notably 
reduce the scatter (e.g. Calzetti et al., 2007; Zhu et al., 2008; Kennicutt et al., 2009). However, 
factors of the combination are different in each work. 

The recently released AKARI/IRC Point Source Catalogue Version /3-1 (hereafter 
IRCPSC) provides positions and fluxes of all-sky survey at S9W (9 /xm) and L18W (18 /xm) 
bands (Ishihara et al., 2010). For numerou MIR data, it would be useful if there is a bench- 
mark of SFR measurement. Comparing with Spitzer 8 and 24 /im bands, the AKARI IRC 
S9W dJid Li 5 bands cover wider wavelength ranges, including silicate absorption features at 
both bands and emission contributed by large PAH molecules at the LlSWhsjid. This work is 
dedicated to investigate whether and to what degree AKARI broadband MIR data could trace 
SFRs, and how much the inclusion of silicate absorption and longer wavelength PAH features 
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would be. 

Start from a sample with multi-wavelength observation flux, we derive the SFR for each 
galaxy by the spectral energy distribution (SED) fitting. MIR data were correlated with the 
SFRs to built the SFR calibrations. Then the results were compared with those from the Spitzer 
observation. The present paper is organized as follows: Section 2 introduces the GALEX-SDSS- 
2MASS-AKARI sample. Section 3 gives a short introduction to methods for calculating SFRs. 
The MIR-SFR relations and some comparisons between AKARI and Spitzer are reported and 
discussed in Section 4. Conclusions are given in Sections 5. 

2. Data 

2.1. Construction of the multi-wavelength sample 

We cross-identified^ IRCPSC with multi-wavelength data constructed by Takeuchi et 
al. (2010) to obtain a sample including MIR photometric measurements. The original sample 
was based on a selection by IRAS-PSCz (Saunders et al., 2000) and AKARI/FIS Bright Source 
Catalog (hereafter FISBSC) data (Yamamura et al., 2008, 2009), which means all these galaxies 
have considerable fluxes at FIR. Then these galaxies were cross-matched with observations by 
GALEX, 2MASS and SDSS. Fluxes at the UV band were obtained by performing specific 
aperture photometry so as not to shred GALEX images. A detailed description of the original 
sample can be found in Takeuchi et al. (2010). 

IRCPSC presents flux data at 9 fim and 18 /im that have an effective bandwidth of 
4.10 yum and a detection limit of 50 mJy, and 9.97 /im and 90 mJy, respectively (Ishihara et al., 
2010). After cross-identification, there were 162 galaxies with a flux of either 9 /im or 18 //m. 
However, nine galaxies were found to have inconsistent fluxes at GALEX, SDSS or 2MASS 
band according to their SEDs, possibly due to the measurement errors or the misidentification 
of objects caused by the inhomogeneous resolution of each observation; they were omitted 
from our sample. We also discarded one galaxy with a too-small redshift {z ~ 0.0008). The 
summary of the sample is listed in Table 1. Searching in the SIMBAD database, most of 
galaxies in our sample are normal star forming galaxies: Only 15 galaxies are classified as 
AGNs (including Seyfert 1 and Seyfert 2 galaxies). All the galaxies in our sample are nearby 
ones. The distribution of the redshift is shown in Figure 1. 

2.2. The Re-estimation of AKARI FIR Data 

The FISBSC flux density for extended sources would be no longer accurate because of 
the point source extraction procedure. In our sample, a considerable fraction of galaxies is 
the extended source. Therefore, it is questionable whether the catalog data of these sources 
are reliable. This is confirmed by making a comparison between FISBSC fluxes and IRAS co- 

^ The searching radius is set at 36" considering the original sample's position accuracy. 
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Fig. 1. The redshift distribution of our sample. The filled area denotes AGNs. 

added fluxes (Figure 2), which were specially calculated for extended sources (Saunders et al., 
2000). In order to investigate the dependence of FISBSC flux on the extension of the galaxy, 
the sample was divided into three sub-samples according to the size of the galaxy (Table 2). 
The lines in Figure 2 give predictions of the difference between two different bands by Dale 
& Helou (2002) (hereafter DH) one-parameter constrained SED templates. The parameter a 
is related to IRAS flux ratio fm/f wo- The log(/6o//ioo) in our sample ranges from —0.60 to 
0.25, and most of the galaxies have log(/6o//ioo) between —0.5 and 0.0, corresponding to an 
a value between 2.625 and 1.375. Accordingly, the model prediction with a = 1.375, 2.625 
and the median value 2.0 are illustrated in Figure 2. Table 3 presents the median values of 
log(/AKARi//iRAs) for each sub-sample. Compared with the model predicted values, there are 
clear discrepancies between FISBSC and IRAS data, especially for 60 — 65 fim and 90 — 100 fim. 
Also note that the large galaxies show greater discrepancies than small ones, indicating the PSF 
photometry is less reliable for larger galaxies, because a considerable part of their flux is left 
out by the relatively small beam size. 

Aiming to obtain more reliable flux data, the photometry of the diffuse maps provided 
by the AKARI group was conducted using SExtractor. The obtained "AUTO" fluxes were also 
compared with IRAS co-added fluxes. 

In the 90 /im band, the result shows a great improvement on the consistency with the 
IRAS flux according to the DH model (Figure 3 and Table 3) for all galaxies in our sample. 
At 65 and 140 /im bands, the consistency is improved for those galaxies in 'Medium' and 
'Large' subsamples, whereas the dispersion increases for the galaxies in 'Small' samples. Less 
improvement is found at 140 /xm band, which can be explained by the relatively larger PSF 
FWHM at this band (~ 60", as compared with ~ 39" for the 65 and 90 /xm bands, Kawada et 
al., 2007). At the 160 nm band, neither fluxes derived from the diffuse maps nor those from the 
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Fig. 2. Comparison between FIS catalog flux and the IRAS flux. Galaxies in different sub-sample are 
shown in different sizes (large open circles: large; small open circles: medium; dots: small). The DH 
models with three different a are shown in different lines (dashed:Q; = 1.375; solid: a = 2.0; dash-dotted: 
a = 2.625). 

IRCPSC are satisfying, which is due to the poor quahty of this band, thus 160 /xm data were 
not used in the following work. Therefore, we kept IRCPSC flux values for 'Small' sources at 
the 60 and 140 /im bands, and applied fluxes derived from diffuse maps to the other sources. In 
addition, fluxes derived from AKARI diffuse maps showed notably smaller measurement errors 
than the IRAS fluxes, greatly improving the data quality. Considering taht the IRAS fluxes 
provide information similar to that given by the AKARI bands, we omitted IRAS fluxes when 
fitting the SEDs. 

Note that the MIR data provided by IRCPSC were little affected by an extension of the 
source, since it applies "AUTO" fluxes by SExtractor, which are suitable for both point sources 
and extended sources. 
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Fig. 3. Comparison between FIS diffuse map flux and tlic IRAS flux. Symbols are of the same meaning 
as in Figure 2 

3. SFR Calculation 

The SED fitting program CIGALE (Noll et al., 2009) was used to calculate the SFR for 
our sample. CIGALE was developed to derive highly reliable galaxy properties by fitting the 
UV/optical SEDs and the related dust emission at the same time, i.e., the stellar population 
synthesized models are connected with infrared templates by the balance of the energy of dust 
emission and absorption. A detailed description of CIGALE can be found in Noll et al. (2009); 
Buat et al. (2010); Giovannoli et al. (2011). Here we give a brief introduction to the main 
features. 

CIGALE allows one to use the stellar SEDs from models given either by Maraston 
(2005) or by Fioc & Rocca-Volmerange (1997, PEGASE). The difference between these mod- 
els is the contribution of the thermally pulsating asymptotic giant branch (TP-AGB) stars. 
In PEGASE models, the contribution from TP-AGB stars is low (Maraston et al., 2006). 
Maraston (2005) increased the contribution from TP-AGB stars adopting the 'fuel consump- 
tion' approach. Maraston et al. (2006) shows that the insufficient consideration of TP-AGB 
stars overestimates the stellar mass by 0.2 dex and worsens the consistency with IR observation 
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data at a redshift of z ~ 2. At lower redshift, it is found that using different models little 
affects the results (Rettura et al., 2006; Eminian et al., 2008). Therefore, the models from 
Maraston (2005) were preferred in this work. The Kroupa IMF (Kroupa, 2001) was used to 
calculate the complex stellar populations (CSPs). 

This code provides two scenarios of star formation: One is "box models" with a constant 
SFR, the other is "r models" , for which SFR decreases exponentially with a typical decay time, 
r. SFR is calculated as SFRbox = Mgai/t for "box models" and SFR^ = Mgai/[r(e*/^ - 1)] for 
"r models" , where Mgai is the galaxy mass (Noll et al., 2009). CIGALE also allows one to apply 
different scenarios for young and old populations. The input SFH here is a constant burst SFH 
for young stellar populations, and an exponentially decreasing one for old stellar populations. 
Thus SFRs were calculated using the formula /ysp ■ SFRbox + /oSP • SFR,-, where /ysp and 
/osp are fractions of young stellar populations and old stellar populations, respectively. 

The attenuation curve adopted by CIGALE is based on a law given by Calzetti et al. 
(2000), with modification of the slope and/or adding a UV bump. The modification of the slope 
is controlled by the factor {X/XvY, i.e., by changing S, the slope of the attenuation curve can 
be modified. We only considered a modification of the slope of the attenuation law here, and 
no bump was introduced. CIGALE allows one to consider different effects of attenuation for 
old and young stellar populations by adding the reduction factor fatt of the dust attenuation for 
the old stellar populations as an input parameter. For the IR part, CIGALE uses DH models, 
which is described in Section 2. Then the dust emission was calculated, and by balancing the 
energy emitted and absorbed, the short and long wavelength parts of the model were connected. 

To compute the output parameters the code provides two methods: "sum" and "max". 
The former calculates the probability distribution functions (PDFs) by taking sums of the 
probability of the models in given bins of parameter space, which might cause an unintentional 
bias when the input parameter values are badly chosen. The latter introduces a fixed number of 
equally sized bins for each parameter and searches the maximum probability of models in each 
bin. Then these maximum probabilities are taken as weights for individual bins to calculate 
the expectation value of the parameter. The advantage of this method is that it alleviates the 
dependence on the choice of parameters (see Noll et al., 2009, for detail). Therefore, we applied 
the "max" method in this work. 

The input parameters applied in this study were adopted from Buat et al. (2010) in 
order to obtain a stable and reliable output. Since CIGALE is unable to trace the unobscured 
emission of an AGN, for Seyfert 1 galaxies, the output decreases the reliability (Buat et al., 
2010). For dust-obscured AGNs, CIGALE provides models to fit the SED, and can avoid 
introducing any severe bias. Therefore, five Seyfert 1 galaxies were rejected in our sample, 
and other AGNs are marked during analysis. In order to keep the accuracy of derived SFRs, 
galaxies with relatively large discrepancies between observation and output spectra (reduced 
> 10) were discarded as being unreliable. At first running, there were 25 such galaxies. By 
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Fig. 4. The comparison between SFRs derived from CIGALE code and from tlie mock galaxies. 

updating the SDSS data using Navigator of SDSS DR7/8 instead of the original pipehne data 
from SDSS DR7, the reduced values of 17 galaxies decreased to less than 10. ^ Four of the 
eight remaining galaxies had very poor quality of AKARI diffuse maps. The other four had 
large extension and brightness, which could cause incomplete flux derivation or saturation in 
the optical plates. Considering that the number was small (~ 6%), these eight galaxies and five 
out of the 17 galaxies were discarded. Therefore, reliable SFRs were derived for 140 galaxies, 
out of which there are 112 with available 9 /im fluxes and 97 with 18 /xm fluxes. 

3.1. The reliability oj the results 

For a SED fitting, the accuracy of the output depends on the input parameters. To give 
proper estimates of the SFRs, the robustness of the results must be tested. A straightforward 
way to check the reliability of the output of CIGALE is to use a sample of mock galaxies with 
comprehensively known physical parameters. Mock galaxies can be generated following a recipe 
in Giovannoli et al. (2011): 1, Run CIGALE on the data of real galaxies. For each galaxies, a 
best model is produced by minimization. Then, from these models, the fluxes at each bands 
can be estimated. 2, Add to each flux a random relative error, which is normally distributed 
with a = 0.1. Thus, we obtain a mock catalog with flux information at every photometric 
band used in this study. The last step is to run the code on the mock catalog and then to 
compare the output parameters with the exact values provided by the best models. The result 
of the comparison is shown in Figure 4. We only present the result concerning SFRs because 
in this work the SFR is the only parameter that needs to be of concern. Figure 4 shows the 
two quantities are well related, indicating SFRs derived here are reliable. 

As discussed in Noll et al. (2009), CIGALE could provide stable results of SFRs as long 
as one constraint beyond PAH band is given. We ran CIGALE with and without MIR data 
to examine the influence on SFR by MIR photometric data. The result shows adding the MIR 

^ There is no modification for other sources which are fitted well. 
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Fig. 5. The comparison between SFR derived witli and witliout MIR data. AGNs are sliown as triangles. 

data or not adding has almost no influence on the resulting value of SFR (Figure 5), while not 
using FISBSC data brings large uncertainties to the output, consistent with the conclusion of 
Noll et al. (2009) using SINGS sample. 

4. Result and Discussion 

Apparently, the luminosity at 9 fim and 18 fim Lg and Lig ^ both correlate with SFRs 

(Figure 6 and Figure 7); the Spearman's correlation coefficients are 0.943 for the Lg-SFR 

relation and 0.956 for the Lig-SFR relation^. Linear regressions following the method provided 

by Kelly (2007) give: 

SFR T 
log -—- = (0.99 ±0.03) log -^-(9.02 ±0.32) (1) 

Mq/jT Lq 

and 

SFR L 

log = (0.90 ±0.03) log —^-(8.03 ±0.30). (2) 

Mq/JT Lq 

The scatters of the data points about the regression lines of Equations 1 and 2 are approximately 
the same, with a = 0.18 dex for 9 fim and 0.20 dex for 18 fim. These tight correlations also 
hold for the surface densities of luminosities and SFRs (Figures 8 and 9). The correlation 
coefficients are 0.961 for 9 fim and 0.945 for 18 fim. The areas of galaxies were calculated from 
5^band images of SDSS. The regression gives: 

log ^ o = (1.02 ±0.03) log ^ 

^ M0yr-%pc-2 ^ ' ^ L^kpc-^ 
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In this paper, L\ refers to vL^, at wavelength A. 



^ Pearson's correlation coefficient is 0.944 for Lg-SFR relation and 0.951 for Lig-SFR relation 
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Fig. 6. The 9 /im luminosity-SFR relation. Tiie daslied line shows the fitting result. The triangles are 
AGNs. 
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Fig. 7. The 18 /im luminosity-SFR relation. Lines and symbols share the same meaning with Figure 6. 

- (9.30 ±0.19) (3) 

with (J = 0.18, and 

(4) 

(5) 
with 0.22. 



log ^^f^ = (0.98 ±0.04) log 

M0yr ^kpc LQkpc 

- (8.89 ±0.25). 



Figures 6 and 7 show that AGNs share a MIR-SFR relation similar to normal star- 
forming galaxies, although the radiation mechanism of AGNs is different from normal galaxies. 
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Fig. 8. The surface densities of 9 /ini luminosity-SFR relation. The dashed line shows the fitting result. 
The triangles are AGNs. 
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Fig. 9. The surface densities of 18 /um luminosity-SFR relation. Lines and symbols share the same 
meaning with Figure 6. 
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Fig. 10. The filter response curves of AKARI 9 fim and 18 /iin bands (dasiied line) and Spitzer 8 /im 
and 24 /im bands (dotted line) . The solid line is the luminosity weighted average spectrum of star forming 
galaxies from Smith et al. (2007). 

The slopes of the logLg-logSFR and logLig-logSFR relations derived here are almost equal to 
one, indicating the MIR-SFR relations are close to linear. 

4.I. Comparison with SFR calibrations from Spitzer data 

Spitzer 8 fim and 24 fim data were investigated as SFR tracers by several authors (e.g 
Wu et al., 2005; Perez- Gonzalez et al., 2006; Calzetti et al., 2007; Relano et al., 2007; Zhu et al., 
2008; Rieke et al., 2009). The 24 /im fluxes were found to be tightly related to the emission of 
the warm dust, and thus more intensively investigated, whereas the 8 /im-SFR relation is more 
complicated, strongly depending on such as the metallicity, size and star-formation history 
(Calzetti et al., 2007), therefore, there are fewer calibrations. 

The main difference between the AKARI and Spitzer filters is their bandwidths. Due 
to the wider band, AKARI 9 /im and 18 /im fluxes are more affected by silicate absorption, 
PAH and molecular hydrogen line emissions (Figure 10). In order to check the reliability of 
our calibrations, the SFRs derived from Equations 1 and 2 were compared with those given by 
Spitzer calibrations. The work to compare was chosen to keep the luminosity range close to the 
present sample (Table 4). The calibrations given by Wu et al. (2005) and Zhu et al. (2008) are 
based on an equation given by Kennicutt (1998) in which Salpeter IMF was used. The use of 
Salpeter IMF will cause smaller SFRs by ~ 0.18 dex than some other IMF with a more shallow 
slope at low masses (Rieke et al., 2009). Therefore this effect was corrected for the results of 
Wu et al. (2005) and Zhu et al. (2008). 

4-1.1. Comparison between SFRs derived from 9 /im and from 8 /im 

Spitzer 8 /xm fluxes were computed from the output spectra of the CIGALE and filter 
response curves of Spitzer IRAC and MIPS. Since the 8 /xm fluxes used in the calibration given 
in Table 4 were the dust emission with the stellar contribution subtracted following the recipe 
of Helou et al. (2004), the 3.6 //m flux was also calculated to compute the stellar composition 
contained in 8 /im, and thus 8 /im dust emission could be obtained (hereafter, we refer to 8 /im 
dust emission as 8 /im emission for conciseness). Note that here the stellar contribution is very 
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Fig. 11. Comparison between tlie SFRs derived from 9 /im emission (Equation 1) and from 8 /im emission 
by Wu et al. (2005) (crosses) and Zhu et al. (2008) (circles). 

small, which can only affect the result by ~ 0.02 dex. The obtained ffux is then converted 
to SFR by formula given by previous work (Table 4). The results are plotted in Figure 11 
(The typical Icr uncertainty for the galaxy of median luminosity is ~ 0.5 dex). The statistical 
information of the comparison is given in Table 5. 

The discrepancy between our results and Wu et al. (2005) may due to several reasons. Wu 
et al. (2005) listed factors such as the accuracy of fiber aperture corrections, the validity of the 
estimation of the obscuration in galaxies by using Balmer decrement, the possible contamination 
to radio and MIR emission from obscured weak AGNs. The larger capacity of our sample (79 
for 8 /im in Wu et al. (2005) compared with 112 in our sample) and the wider coverage of 9 /xm 
band may also cause such difference. Another possible reason is that the oversimplification of 
PAH emission in DH models underestimates the 8 yum flux and therefore gives smaller SFRs. 
However, this level of discrepancy is well within the scatters in Equations 1. 

The discrepancy between Wu et al. (2005) and Zhu et al. (2008) is because Zhu et al. 
(2008) included 8 /im-weak Hll galaxies with lower MIR luminosity (Zhu et al., 2008). Since no 
such galaxies were included in present sample, it is reasonable that present result agrees with 
Wu et al. (2005) better. 

4-1-2. Comparison between SFRs derived from 18 /im and 24 /im 

For a 18 yum-24 yum SFRs comparison, more reference calibrations are available (Table 
4). The 24 /im fluxes are also derived from the output spectra of the SED fitting by CIGALE. 
The statistical information of the comparison is given in Table 6. The converted SFRs are 
plotted in Figure 12. Since the SFRs derived from Wu et al. (2005) are quite similar to those 
from Zhu et al. (2008) and the luminosity range in Wu et al. (2005) is closer to this work, the 
results from Zhu et al. (2008) are omitted for the sake of brevity. 
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Fig. 12. Comparison between tlie SFRs derived from 18 ^m (Equation 2) and from 24 fim emission by 
Wu et al. (2005) (crosses), Calzetti et al. (2007) (triangles) and Rieke et al. (2009) (dots). The dotted line 
gives the lower limit where the calibration of Rieke et al. (2009) applies. 

The results of Wu et al. (2005) and Zhu et al. (2008) agree well with our result after 
the correction of IMF. Rieke et al. (2009) assembled SED templates for local luminous and 
ultraluminous infrared galaxies and combined the result of Dale et al. (2007) and Smith et al. 
(2007) to produce templates at lower luminosities. Their result is applicable to galaxies with 
24 /xm luminosity greater than 6 x 10* Lq, corresponding to log SFR = —0.33, which is 
shown by the dotted line in Figure 12. The present result agrees very well with Rieke et al. 
(2009) above the hmit. 

Our result is a little higher than the one given by Calzetti et al. (2007). A possible 
reason is that the result of Calzetti et al. (2007) was derived for HII clouds by Paa emission, 
which might be poorly applied to galaxy-wide calculations, because the diffuse MIR or Pao; 
emissions in the whole galaxy are not included (Alonso-Herrero et al., 2006; Calzetti et al., 
2007; Kennicutt et al, 2007; Rieke et al., 2009). 

4-2. Combination of FUV and MIR indicators 

At lower IR luminosity, the IR indicator may fail to trace part of the UV photons 
from young stars due to the increased transparency of the ISM. A combination of unobscured 
FUV and MIR luminosities, (FUV + aMIR), may efficiently compensate for the lost energy, 
and could trace the SFR linearly (Zhu et al., 2008). However, upon converting SFRs to dust 
obscuration corrected FUV fluxes by eq.[l] from Kennicutt (1998), we find that a = 2.53 and 
3.33 with a scatter of 0.17 dex for 9 fim and 0.20 dex for 18 /im, respectively. The scatter is not 
reduced significantly. This fact indicates that the origin of the scatter in MIR-SFR diagram 
is complicated: not only the untraced UV photons, but also other unknown factors, such as 
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the variation of the physical conditions within each galaxy, the distribution of dust and photo 
dissociation regions (PDRs), etc. 

4.3. Metallicity 

We attempt to investigate the gas-phase metallicity range of our sample by searching in 
the metallicity database measured by Tremonti et al. (2004) for SDSS galaxies. Unfortunately, 
the 12 + log(0/H) values are given for only 33 galaxies (all higher than 8.75). Therefore, we 
applied a compromised method: to investigate the stellar mass of our sample. Measured 
by CIGALE, the stellar mass M^, of all the galaxy in the sample is larger than 10®'^ Mq. Thus, 
from mass- metallicity relation given by Tremonti et al. (2004): 



the metallicity 12 + log(0/H) for our sample is higher than 8.4. This is not surprising because 
the initial sample from Takeuchi et al. (2010) is IR selected, which means a considerably high 
luminosity in IR and thus sufficient dust content and relatively high metallicity range. Thus, the 
MIR-SFR relations derived here could only be extrapolated to other high-met allicity galaxies. 
The situation for low metallicity galaxies is rather complicated; since the opacity of the galaxy 
decreases, MIR would be unable to trace most of the UV/optical photons and thus lose the 
ability as an SFR indicator (Calzetti et al., 2007). 

4.4. AGNs 

Although AGNs have distinct features from normal galaxies in various physical prop- 
erties, they share the same trend in the MIR-SFR figures. A possible reason is that the 
contribution from AGN component is minor (from SED fitting, less than 15%), therefore the 
host galaxy component dominates the spectrum. To investigate the effect of AGNs on MIR 
emission, we plot average SEDs for AGNs and normal galaxies in Figure 13. On average, AGNs 
are brighter than normal galaxies at all bands, while with a lower MIR/90 /im ratio (Table 7). 
There are two possible reasons: 1, silicate absorption occurs more strongly in AGNs. 2, PAH 
molecules are destroyed by the harsh radiation field in AGNs. The latter is more convincing. 
Studies show that small PAH molecules contributed to shorter wavelength MIR emission are 
destroyed more easily than large ones contributed to longer wavelength MIR emission (Smith 
et al., 2007; Treyer et al., 2010, and reference therein), which is consistent with the lower 9 /im 
flux value than 18 /xm in Figure 13. Nevertheless, these differences between AGNs and normal 
galaxies are only on their average level, which could not be used to distinguish AGNs and 
normal galaxies in individual cases. 



12 + log(0/H) 



1.492 + 1.847(logM,) 
0.08026(logM,)2, 



(6) 



15 



2 r 1 I 1 !| p !l 




-4 E ii '■■ ii 

0.1 1.0 10.0 100.0 



Fig. 13. The average SED (dashed hnes) and the SED normaUzed at 90 /im (sohd hnes) of AGNs (sohd 
symbols) and normal galaxies (open symbols). 

5. Conclusion 

We combined AKARI/IRC 9 yum and 18 /xm data with a previous sample to construct a 
FIR selected multi-wavelength sample with MIR photometric measurements. The FIR data of 
AKARI/FIS in the original sample were re-estimated by photometry of AKARI diffuse maps 
to correct the bias of PSF photometry for extended sources. Then, the SEDs of the sample 
were fitted by CIGALE, and the SFRs were obtained. Regression analysis was conducted to 
investigate MIR-SFR relations. SFRs converted from AKARI MIR fluxes were compared with 
those from the Spitzer MIR fluxes to test the reliability of AKARI MIR-SFR calibrations. From 
the result, we draw the following conclusions: 

1. Both 9 /im and 18 /im luminosities correlate with SFRs, and thus could be converted 
SFRs. 

2. A combination of FUV and MIR luminosity barely reduces the scatters, indicating that 
the unobscured UV photons are not the only reason of the variation of MIR-SFR relation. 

3. A comparison of the SFRs derived from Equations 1 and 2 with the ones derived from 
Spitzer MIR-SFR relations shows that the silicate absorption included in S9W {9 fim) and 
L18W (18 /im) bands little affects the results. The discrepancies, if any, are well within 
the uncertainties. 

4. AGNs in the sample show no discrepancy with normal galaxies in the MIR-SFR diagrams. 
The smaller average MIR fluxes for AGNs than normal galaxies might indicate the small 
PAH molecules are destructed by harsh radiation from AGNs. 

In summary, for IR selected galaxies the rest frame 9 /im and 18 /im emissions are efficient 
tracers of SFRs, and the equations derived here should be applicable to other dust rich galaxies. 
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Table 1. A brief summary of the sample. Data with an asterisk (*) are not used for SED fitting in Section 3. 



Survey Band Wavelength (^m) Nb. of sources 

GALEX FUV, NUV 0.153, 0.231 153 

SDSS u, <?, r, i, z 0.355, 0.469, 0.617, 0.748, 0.893 153 

2MASS J, E, Ks 1.244, 1.655, 2.169 153 

IRAS* band-1, 2, 3, 4 12, 25, 60, 100 153 

AKARI IRC S9W 9 126 

AKARIIRC L18W 18 106 

AKARI FIS N60, Wide-S, Wide-L, N160* 65, 90, 140, 160 153 



Table 2. Criteria to divide the sample into three sub-samples. The length of major axis a and minor axis b of each galaxy 
is obtained from the SDSS image. The 40" is taken as the separating value considering the PSF size of AKARI/FIS. 

Sub-sample Criteria 

Small a < 40" and b < 40" 

Medium a > 40" or 5 > 40" 

Large a > 40" and b > 40" 



Table 3. Comparison between AKARI and IRAS flux. The mean values of the difference between AKARI and IRAS fluxes 
according to AKARI FISBSC fluxes (Catalog), fluxes derived from diffuse maps (Map) and DH model prediction (Model) 
are listed for each subsample of galaxies ('L', 'M' and 'S' represent large, medium and small group in Table 2, respectively). 







Catalog 






Map 






Model 






L 


M 


S 


L 


M 


s 


a = 1.375 


a = 2 


a = 2.625 


log(/65/./60) 


-0.18 


-0.09 


-0.03 


-0.10 


-0.12 


-0.01 


0.00 


0.03 


0.04 


log(/9o//l0o) 


-0.33 


-0.27 


-0.24 


-0.13 


-0.12 


-0.11 


0.05 


0.00 


-0.04 


Iog(/l4o//l0o) 


-0.09 


-0.04 


-0.09 


0.03 


0.01 


-0.02 


-0.30 


-0.11 


0.02 


Iog(/l6o//l0o) 


-0.20 


-0.17 


-0.32 


0.01 


0.00 


0.07 


-0.45 


-0.20 


-0.02 
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Table 4. SFR calculations based on Spitzcr data. The luminosity is expressed in Lq and SFR in A/Q/yr. 



Work 


SFR calculation 


Wu et al. (2005) 


log SFR 


= (1.09 ± 0.06)logL8 - (10.03 ± 0.16) 


Zhu et al. (2008) 


log SFR 


= (0.93 ± 0.03)logL8 - (8.59 ± 0.08) 


Wu et al. (2005) 


log SFR 


= (0.89 ± 0.06)logL24 - (7.82 ± 0.17) 


Calzetti et al. (2007)* 


log SFR 


= 0.88501ogL24 - 8.17 (1 a uncertainty 0.03) 


Zhu et al. (2008) 


log SFR 


(0.85 ± 0.01)logi24 - (7.47 ± 0.06) 


Rieke et al. (2009) 


SFR = 
SFR 


7.8 X 10-i"i24 for 6 X 10^ < L24 < 1-3 x 10^° Lq 

7.8 X 10-i"i24 X (7.76 X 10-"i24)° °'^^ for ^24 > 1.3 x 10^" L© 


* For galaxies with 12 + 


log{0/H) 


> 8.35, i.e., "high mctallicity" galaxies in Calzetti et al. (2007). 



Table 5. Regressions between SFRs derived from AKARI 9 /im and from Spitzer 8 jrar. 
logSFR(8 /xm) = a + 61ogSFR(9 tira) and the mean value of the difference: < logSFR(8 puxa) -logSFR(9 /^m) >. 



8 /im calibration 


a b 


Correlation coefficient 


<logSFR(8 fim) 


-logSFR(9 Mm) > 


Wu et al. (2005) 


-0.17 ± 0.02 1.05 ± 0.03 


0.968 


-0.14 


± 0.18 


Zhu et al. (2008) 


-0.20 ± 0.02 0.90 ± 0.02 


0.968 


-0.25 


± 0.16 



Table 6. Regressions between SFRs derived from AKARI 18 fim and from Spitzcr 24 fj,m: 
logSFR(24 /xm) = a + 61ogSFR(18 fim) and the mean value of the difference: < logSFR(24 /^m) -logSFR(18 /xm) >. 



24 calibration 




a 






b 




Correlation coefficient 


< logSFR(24 /ini) 


-logSFR(18 /im) > 


Wu et al. (2005) 


-0.02 


± 


0.02 


1.00 


± 


0.02 


0.982 


-0.01 


± 0.12 


Calzetti et al. (2007) 


-0.23 


± 


0.02 


1.00 


± 


0.02 


0.982 


-0.23 


± 0.12 


Zhu et al. (2008) 


-0.02 


± 


0.02 


0.96 


± 


0.02 


0.982 


-0.05 


± 0.11 


Rieke et al. (2009) 


-0.14 


± 


0.02 


1.14 


± 


0.02 


0.982 


-0.04 


± 0.17 



Table 7. The MIR to 90 fim flux ratio for AGNs and normal galaxies. 



Gal. 


log(/9//9o) 


log(/l8//90) 


l0g(/9//l8) 


AGN 


-1.69 


-1.35 


-0.34 


Normal 


-1.35 


-1.19 


-0.16 
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